Heat transfer in lattice BGK modeled fluid

The thermal lattice BGK model is a recently suggested numerical tool aiming at solving problems of thermohydrodynamics. The quality of the lattice BGK simulation is checked in this paper by calculating temperature profiles in the Couette flow under different Eckert and Mach numbers. A revised lower order model is proposed to improve the accuracy and the higher order model is proved to be advantageous in this respect, especially in the flow regime with a higher Mach number.     

Multibody dynamics of aircraft occupants seated behind interior walls

To improve aircraft crash safety, conditions critical to occupants survival during a crash must be known. In view of the importance of this problem, studies of post-crash dynamic behavior of victims are necessary in order to reduce severe injuries. In this study, crash dynamics program SOM-LA/TA, incorporating a dynamic model of the human body with a finite element model of the seat structure was used. Modifications were performed in the program for reconstruction of an occupant's head impact with the interior walls or bulkhead. A viscoelastic-type contact force model was used to represent the compliance characteristics of the bulkhead. Correlated studies of analytical simulations with impact sled test results were accomplished. A parametric study of the coefficients in the contact force model was then performed in order to obtain the correlations between the coefficients and the Head Injury Criteria. A measure of optimal values for the bulkhead compliance and displacement requirements was thus achieved in order to keep the possibility of a head injury as little as possible. This information could in turn be used in the selection of suitable materials for the bulkhead, instrument panel, or interior walls of an aircraft.     

Molecular dynamics simulation of entropy and surface tension for grain boundary of α-Fe

The grain boundary is an interface and the surface tension is one of its important thermodynamic properties. In this paper, the surface tension of the ∑9 grain boundary for α-Fe at various temperatures and pressures is calculated by means of Computer Molecular Dynamics (CMD). The results agree satisfactorily with the experimental data. It is shown that the contribution of entropy to surface tension of grain boundary can be ignored. The project supported by the National Natural Science Foundation of China and the Science Foundation of Chinese Academy of Sciences.     

A numerical study of time-dependent schr?dinger equation for multiphoton vibrational interaction of NO molecule, modelled as Morse oscillator, with intense far-infrared femtosecond lasers

For the NO molecule, modelled as a Morse oscillator, time-dependent (TD) nuclear Schr?dinger equation has been numerically solved for the multiphoton vibrational dynamics of the molecule under a far-infrared laser of wavelength 10503 nm, and four different intensities,I = 1 × 10⁸, 1 × 10¹³, 5 × 10¹⁶, and 5 × 10¹⁸ W cm⁻² respectively. Starting from the vibrational ground state at zero time, various TD quantities such as the norm, dissociation probability, potential energy curve and dipole moment are examined. Rich high-harmonics generation (HHG) spectra and above-threshold dissociation (ATD) spectra, due to the multiphoton interaction of vibrational motions with the laser field, and consequent elevation to the vibrational continuum, have been obtained and analysed. Dedicated to Professor C N R Rao on his 70th birthday An erratum to this article is available at .     

Robust Dynamics and Control of a Partially Observed Markov Chain

In a seminal paper, Martin Clark (Communications Systems and Random Process Theory, Darlington, 1977, pp. 721–734, 1978) showed how the filtered dynamics giving the optimal estimate of a Markov chain observed in Gaussian noise can be expressed using an ordinary differential equation. These results offer substantial benefits in filtering and in control, often simplifying the analysis and an in some settings providing numerical benefits, see, for example Malcolm et al. (J. Appl. Math. Stoch. Anal., 2007, to appear). Clark’s method uses a gauge transformation and, in effect, solves the Wonham-Zakai equation using variation of constants. In this article, we consider the optimal control of a partially observed Markov chain. This problem is discussed in Elliott et al. (Hidden Markov Models Estimation and Control, Applications of Mathematics Series, vol. 29, 1995). The innovation in our results is that the robust dynamics of Clark are used to compute forward in time dynamics for a simplified adjoint process. A stochastic minimum principle is established.     

On the Use of the Restitution Condition in Flexible Body Dynamics

The difference between the classical treatment offlexible body impact and the treatment of impact in flexiblemultibody dynamics is due to several fundamental reasons. Inthe classical impact theory, simple structures such as beamsand plates are used. Infinite dimensional models can bedeveloped for these simple structural elements to study theimpact dynamics and the wave propagation problem. Flexiblemultibody impact problems, on the other hand, involve bodieswith complex geometry that cannot be modeled using infinitenumber of degrees of freedom. Furthermore, the classicalimpact theory has been mainly concerned with the impactbetween a rigid mass that moves without constraints beforeit impacts a simple flexible structure. This is not amultibody simulation scenario in which the impact occursbetween kinematically constrained bodies that are subjectedto impulsive constraint forces in addition to the impactforces. These constraint forces can influence the motion ofthe two bodies immediately after impact, and as aconsequence, the simple classical theory scenario of impactdoes not apply. It is the objective of this paper to discussthe use of the restitution condition in flexible multibodyimpact problems and demonstrate that the use of thisapproach does not exclude the classical formulation.Nonetheless, the impulse momentum balance approach can serveas an effective and efficient procedure for solving theimpact problem in finite dimensional models that do not obeythe classical wave theory. Energy results of simplestructural elements are presented in order to demonstratethe consistency of using the impulse momentum balanceapproach in solving impact problems in finite dimensionalflexible body applications.     

The effect of nanoparticle shape on polymer‐nanocomposite rheology and tensile strength

Nanoparticles can influence the properties of polymer materials by a variety of mechanisms. With fullerene, carbon nanotube, and clay or graphene sheet nanocomposites in mind, we investigate how particle shape influences the melt shear viscosity η and the tensile strength τ, which we determine via molecular dynamics simulations. Our simulations of compact (icosahedral), tube or rod‐like, and sheet‐like model nanoparticles, all at a volume fraction ? ≈ 0.05, indicate an order of magnitude increase in the viscosity η relative to the pure melt. This finding evidently can not be explained by continuum hydrodynamics and we provide evidence that the η increase in our model nanocomposites has its origin in chain bridging between the nanoparticles. We find that this increase is the largest for the rod‐like nanoparticles and least for the sheet‐like nanoparticles. Curiously, the enhancements of η and τ exhibit opposite trends with increasing chain length N and with particle shape anisotropy. Evidently, the concept of bridging chains alone cannot account for the increase in τ and we suggest that the deformability or flexibility of the sheet nanoparticles contributes to nanocomposite strength and toughness by reducing the relative value of the Poisson ratio of the composite. The molecular dynamics simulations in the present work focus on the reference case where the modification of the melt structure associated with glass‐formation and entanglement interactions should not be an issue. Since many applications require good particle dispersion, we also focus on the case where the polymer‐particle interactions favor nanoparticle dispersion. Our simulations point to a substantial contribution of nanoparticle shape to both mechanical and processing properties of polymer nanocomposites. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1882–1897, 2007     

Capacity augmentation of a supply chain for a short lifecycle product: A system dynamics framework

For an innovative product characterized by short product lifecycle and high demand uncertainty, investment in capacity buildup has to be done cautiously. Otherwise either the product’s market diffusion is impeded or the manufacturer is left with unutilized capacity. Using the right information for making capacity augmentation decisions is critical in facing this challenge. In this paper, we propose a method for identifying critical information flows using the system dynamics model of a two-echelon supply chain. The fundamental premise of system dynamics methodology is that (system) structure determines (its) behavior. Using loop dominance analysis method we study the feedback loop structure of the supply chain system. The outcome is a set of dominant loops that determine the dynamics of capacity growth. It is revealed that the delivery delay information has little effect while the loop that connects retail sales with production order affects the dynamics significantly. Modifying this loop yields appropriate capacity augmentation decisions resulting in higher performance. What-if analyses bring out effects of modifying other structural elements. In conclusion, we claim that the information feedback based methodology is general enough to be useful in designing decision support systems for capacity augmentation. The limitations of the model are also discussed and possible extensions identified.     

General synchronization dynamics of coupled Van der Pol–Duffing oscillators

This paper considers the general synchronization dynamics of coupled Van der Pol–Duffing oscillators. The linear and nonlinear stability analysis on the synchronization process is derived through the Whittaker method and the Floquet theory in addition to the multiple time scales method. A stability map displaying different dynamical states of the system is performed. Numerical simulation is carried out to support and to complement the accuracy of the analytical treatment.